POLYLOGARITHMS, HYPERFUNCTIONS AND 
GENERALIZED LIPSCHITZ SUMMATION FORMULAE 



STEFANO MARMI AND PIERGIULIO TEMPESTA 



Abstract. A generalization of the classical Lipschitz summation for- 
mula is proposed. It involves new polylogarithmic rational functions 
constructed via the Fourier expansion of certain sequences of Bernoulli- 
type polynomials. Related families of one-dimensional hyperfunctions 
are also constructed. 
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1. Introduction 

The purpose of this paper is to provide a natural setting which allows to 
generalize the Lipschitz summation formula to negative powers. 

We recall that the classical Lipschitz summation formula gives the Fourier 
series expansion of the periodic analytic function obtained by summation 
over integer translates of the power , where z € ^ (the complex upper 
half plane) and A: is a positive integer. 

By introducing suitable generalizations of the Lipschitz formula, we will 
show how to construct new classes of hyperfunctions. Indeed, each of these 
generalized formulae will provide a hyper functional equation involving a 
specific two-variable polylogarithm series. 
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A simple introduction to the theory of hyperfunctions is sketched in Sec- 
tion 2, where the interested reader is referred to the classical works |15j . [7]. 
Hyperfunctional cohomology groups have also been related to automorphic 
forms and period functions [3]. 

The main results of this paper are essentially two. The first result is a 
generalization of the classical Lipschitz summation formula. The second one 
is the connection we establish between a new class of polylogarithms, the 
theory of hyperfunctions in one variable and that of polynomial structures 
of Appell type. 

Precisely, our construction can be summarized as follows. We will define 
suitable sequences of Appell polynomials of Bernoulli type: indeed, they 
share with the Bernoulli polynomials several arithmetic properties, includ- 
ing certain famous congruences (see the Appendix). Their periodic versions 
provide primitives of the periodic delta function. To each sequence we will 
associate a natural extension of the notion of polylogarithmic function de- 
fined in the unit disc. We call this extension a delta rational function. It 
is a two-variable Dirichlet series, extending to the whole Riemann sphere 
as a meromorphic function. The case of the classical Bernoulli polynomials 
corresponds to the standard polylogarithmic function, and it is treated thor- 
oughly. The new Appell polynomial sequences we construct, {Pn (x)}^^^ 
and {Qn (depending on certain parity properties) and the general- 

ized polylogarithms (5„ {q) proposed here satisfy interesting hyperfunctional 
equations. These equations in turn correspond to generalized Lipschitz sum- 
mation formulae. For instance, in the case of polynomials {P„ (3^)}ngN it 
reads 

(1) 

f A_„ (q) if \q\ < 1, i.e. > 

^ifp^ (r + fc) =2i(27ri)-" <^ 

feez " [ (-1)""^ A_„ (g-i) if\q\>l, i.e. < 

where Pn are the hyperfunctions associated to the polynomial P„ and ipp 
is the function in (C\ [0, 1]) representing 

In the last part of the paper, the relation between the theory of formal 
groups and that of hyperfunctions is clarified. We will show that the Appell 
structures introduced in the previous construction can be viewed as poly- 
nomial realizations of certain formal group laws, related with the universal 
Lazard formal group. 

The future research plans include an extension of the proposed construc- 
tion to the multidimensional case (see also |12] for a different generaliza- 
tion), as well as possible applications of the proposed Lipschitz formulae to 
the study of Eisenstein series and periods of modular forms. 

Acknowledgements. One of us (P. T.) is grateful to U. Zannier for a 
very useful discussion. 
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2. Hyperfunctions: basic preliminaries 

In this section, we will provide a brief and self-consistent introduction 
to the theory of hyperfunctions of a single variable, following closely [15j , 
Chapter IX of [7|, [10] and [9]. For a more extensive treatment and fur- 
ther details, the reader is invited to consult these books, as well as [13] for 
interesting applications. 

Let us denote by O the sheaf of holomorphic functions on C, and denote 
by C+ and C~ the upper and lower half planes of C. 

Definition 1. The space of hyperfunctions B on the real line M is 

(2) Bm:=Hi{C,0), 
i.e. is the first sheaf cohomology group on M. 

Now, since U = C \ M, we have the following decomposition: 

(3) Hi (C, O) = [H^ (C+, O) © (C-, O)] / (C, O) . 

In other words, since (O) represents nothing but the global sections of 
the sheaf (i.e. holomorphic functions), a hyperfunction can be thought as 
a pair of holomorphic functions on the upper and lower half planes respec- 
tively, modulo an entire function. This geometric definition immediately 
generalizes to the case of hyperfunctions in several variables. If C M is an 
open set, and U an arbitrary complex neighborhood of C/, then clearly 

(4) Bm = H},{U,0). 

Still denoting by O (U) the space of holomorphic functions in U, another 
equivalent definition is 

(5) B{n) -.= 11111 o{u\n) / oiu) 

where the inductive limit with respect to the family of complex neighbor- 
hoods [/ D is considered. A hyperfunction / (x) is therefore an equivalence 
class [F (z)], whose representative is F (z) O {U \ Q). The representative 
F (z) is said to be a defining function of / (x). Since 

u Dn = u+uU-, c/± = c/n{9z^o}, 

often the following boundary-value representation is used: 

/ (x) = F+ (x + iO) - F_ (x - iO) , 

withF±(z) = F{z)\^^. 

Alternatively, one can construct a theory of hyperfunctions based on an- 
alytic functionals. Let K C C be a non empty compact set, and denote by 
A the space of entire analytic functions in C. 

Definition 2. The space A' (K) of the analytic functionals carried by K is 
the space of linear forms u acting on A such that for every neighborhood V 
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of K there is a constant Cy > sueh that 

(6) \u{ip)\<Cv s\xp\(p\, G A. 

V 

Observe that A' (K) is a Frechet space, since a seminorm is associated to 
each neighborhood V of K. One can define 

(7) B{n):=A'{n)/A'{dn). 

It is interesting to notice that the space of hyperfunctions tp G B (Q) with 
compact support K G ^ can be identified with analytic functionals in A' (M) 
with support K. Indeed, an analytic functional u on U^^^Ki can always be 
decomposed into a sum u = ui + ... + Uj-, with each of the functionals 
Uj G A' (Kj). Consequently, since supp C K U dQ, the contribution of 
ijj on can be factored out and tp is identified with an uniquely defined 
functional with support in 0,. 

Therefore, we can also think of A' (K) as the space of hyperfunctions 
with support in K. The link between the two approaches to the theory 
of hyperfunctions is now provided by the following lemma. Let 0^(C \ 
K) denotes the space of holomorphic functions on C \ K and vanishing at 
infinity. 

Lemma 3. The spaces A' (K) and 0^{C \ K) are canonically isomorphic. 
To each u € A' (K) it corresponds a function ip £ 0^{C \ K) given by 

<f{z)=u{c^), \/zeC\K, 

where Cz (x) = ^^zj- Conversely, to each (p & 0^{C \ K) it corresponds 
the hyperfunction 

(8) u{i^) = ^ j^p>{z)^{z)dz, VV'G^, 

where 7 is any piecewise C^path winding around K in the positive direction. 

For future purposes, we also briefly describe periodic hyperfunctions. Let 
T^ = R\ ZcC\Z. A hyperfunction on is a linear functional ^' on the 
space O (T^) of functions analytic in a complex neighborhood of such 
that for all neighborhood V of there exists a constant Cy > such that 

(9) |^'($)| < sup|$|, V$gO(F). 

V 

We will denote by A' (T^) the Frechet space of hyperfunctions with sup- 
port in T^. Let Os denote the complex vector space of holomorphic functions 
$ : C \ R — > C, 1-periodic, bounded at ± ioo and such that $ (±zoo) := 
limcj2;_^-l-oo ^ (z) exist and verify $ {+ioo) = — $ {—ioo). 

Lemma 4. The spaces A' (T^) and O^, are canonically isomorphic: to each 
^ & A' (T^) it corresponds a function $ G Os given by 

^{z) = ^{Cz), yzeC\K, 
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where Cz {x) = cot tt {x — z). Conversely, to each $ G it corresponds 
the hyperfunction 

(10) ^{~)='-j^^{z)^{z)dz, VHeA'(Ti), 

where V is any piecewise path winding around a closed interval / C K 0/ 
length 1 in the positive direction. 

Given a compactly supported hyperfunction, making infinitely many copies 
of it translating its support leads to a periodic hyperfunction. This point of 
view is systematically exploited in [U], from which the following commutative 
diagram is taken: 

^'([0,1]) Oi(C\[0,l]) 

(11) A'{T^) 

The horizontal lines are the above mentioned isomorphisms, and Ylz 
notes the summation over integer translates. 

3. Delta rational functions 

3.1. Delta rational functions and polylogarithms. In this section we 
introduce the notion of delta rational functions, which represent a natural 
generalization of the notion of polylogarithms [11], [19]. These functions are 
connected to the classical Bernoulli polynomials via the Fourier expansion 
formula (jl4p and enter a generalized Lipschitz summation formula providing 
new classes of hyp erf unctions. 

Definition 5. The delta rational function is the series expressed for any 
n G Z, and q GO by 

00 

(12) <5n(g)= J^feV- 

k=l 

Observe that, if n > 0, 5n extends to the whole Riemann sphere as a 
rational function of degree n + 1 with just a simple pole of order n + 1 at 
q = 1. Ifn<— 1 then 6n coincides with the classical polylogarithm series 
of order —n (actually 6-1 (q) = — log(l — q)) and extends to the whole C 
\[1, +00) and as a multiplicative function to the whole C \{0, 1, 00}. Indeed, 
since 

qdqSn (q) = 6n+i (q) , 
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one can define the analytic continuation of For instance, the continuation 
of 6-2 is obtained by means of the integral formula 



t Jo \h 1-C/ t 

Note that [1, +oo) is a branch cut. For all n E Z and g E D, one has 

(13) 5n{q'')=k-'-'' Yl '^"(^'?)' 

Afc=l 

where A denotes a fe-th root of unity. The equality extends to the closed 
disk if n < —2. One can directly prove the following result. 

Lemma 6. The fundamental inversion equation holds 

'0 if n>l 

(14) Sniq) + i-ir6n{q-') = \ 

for all q ^ 1 if n > 0, and q in C\ [0, +oo] if n < 0. Here is the k-th 
Bernoulli polynomial. 



3.2. Periodic hyperfunctions. The inversion relation (jl4p has a beautiful 
interpretation in terms of hyperfunctions, as we will show below. 

Definition 7. The periodic Bernoulli functions and distributions are ex- 
pressed by 



(15) Bn (x) 



Bn{x-[x]) 



1,! 



if n > 1 



— 1 + 5j {x) if n = , 
4""^ {x) if n<0 
where [x] is the integer part of x, 5j is the periodic delta distribution and 

(k) 

6 J its derivative of order k >0. 
Explicitly, 

+ 00 

(16) St{x)= e^"*'''- 

fc=— oo 

From the well-known property of the periodic Bernoulli functions 

(17) £^{B^{x-[x])) = n\{l-6j), 

we see that we can consider the Bernoulli functions as primitives of the 
periodic delta function. It is immediate to check that the following hyper- 
functional equation holds. 

Proposition 8. Let q^ = ^'^'^^(^^^0) ^ a; E M. For all n ^ Z we have 

(18) 5„(g+) + (-l)"5„((g-)"') =(2^i)-"S_„(x). 
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Note that on both sides one can apply derivatives qdq and (27ri) ^ dx- 

3.3. Generalized Lipschitz summation formula. As is well known, the 
classical Lipschitz formula is a consequence of the Poisson summation for- 
mula, relating sums over integers of pairs of Fourier transforms. It states 
that 



(for a proof of ([19]) . see, for instance, D. Zagier, Chapter 4 in p8|). 

In this section, we will give an hyper functional generalization of ()19p . We 
make use of the fact that periodic hyperfunctions with compact support 
are the result of a summation over integer translates of hyperfunctions with 
support on [0, 1], according to the commutative diagram (jlip . By applying 
the results described in the previous section, we get a functional version of 
the Lipschitz formula. 

Theorem 9. For all n G 7^ we have 
(20) 

6-n (q) if \q\ < 1, i.e. 9r > 



Y,VBjr + k)=2i{27ri)-'' < 



[ (-l)"-^5_„,(g-i) if \q\ > 1, i.e. < 



whereas the usual Lipschitz formula corresponds to n < —1. Here Bn is the 
restriction to [0,1] of Bn and ip-^^ is the function in (C\[0, 1]) which 
represents the hyperfunction Bn 

(21) ^(T) = (Br, ^ ^ 



TT X — T 



[0,1] 



Proof. It is enough to remark that an inverse of the operator on the hy- 
perfunctions is just the restriction of a periodic hyperfunction to the interval 
[0, 1]. Thus the r.h.s. of ^ reads 

(22) j;(27ri)-"B_„(x + A;) 

fcez 

and from (jl4p and (jlSp one obtains the desired formula by considering the 
two associated holomorphic functions. □ 

Explicitly, formula (|2ip reads 
(23) 



[0,1] 

1 + 5 (x) , = -i^[^ + r log (1 - 1/r)] ^fn = 

-^Bn (x) , i^)^^ ^^ = -^. [Bn (r) log (1 - 1/r) + i?„ (r)] ifn>0 
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where Rn-i is the polynomial of degree n — 1 such that (r) log (1 — 1/r) + 
Rn (t) € (C\ [0, 1]). Here are the first six polynomials 

RiiT) = l, R,{t) = t-^, R,{r)=T'-T+^, 

R^ (^) = ^3 _ 3^2 1^ 1 ^ ^ ^4 _ + 3^2 + _ ii, 

^ ^ 2 3 12 ^ ^ 4 4 360 

(24) i?6 (r) = t5 - + + -r^ _ l^r - — . 
^ ^ ^ ^ 2 3 2 60 120 

Remark. One could object that we should have used the simplest Appell 
sequence and their associated hyperfunctions on [0, 1] to write the 

generalized Lipschitz summation formula. But this wold trivially give zero, 
since 

(25) = ^— [Bn+l (x + 1) - Bn+i (x)] 

n + 1 

and the function Bn+i (x + 1) — B^+i (x) clearly belongs to the kernel of 
Y^^- Instead, other less trivial Appell sequences can be used to provide 
useful generalizations of the construction we proposed, as will be illustrated 
in the subsequent sections. 



4. The general case: sequences of Appell polynomials and 

hyperfunctions 

4.1. Appell sequences. By analogy with the theory developed in the pre- 
vious section, we study a more general class of polynomials of Bernoulli 
type, with the aim of constructing associated families of hyperfunctions. 
The main request is that they belong to the class of Appell polynomials, i.e. 
polynomials satisfying the property 

(26) -^An (x) = n An-i (x) , 

ax 

with the normalization 

(27) Aq (x) = const. 

Lemma 10. The Fourier series expansion of the polynomials 126\) . for < 
X < 1 and n > 1 has the form 

oo 

(28) An{x)= ^ c,(n)e2-^^ 

fc=— oo 

with 

(29) Ck = C An{t)e-^^''''dt. 

Jo 
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We get 
(30) 



oo 



An (x) = -n! 



k=l 



'fij 2TTikx , /n^+l-i „-27Tikx , 



wzi/i < a; < 1, and c/^^- = Aj (1) — ^-, (0), j = l,..,n. The standard 
Bernoulli polynomials correspond to the case cq = 0, = 1, <Pj = 0, 
j = 2,3,... 

Proof. It is a direct consequence of the conditions (|26p - (|29|) and of the for- 
mula of integration by parts. □ 

If ipj = for j even or j odd, then formula (jSOp can be written in a more 
compact form. The corresponding polynomial sequences will be denoted by 
{Pn {x)}nm iln (^)}ngN' respectively. 

a) If ipj = for j even, we define 
(31) 



1 fi 



j odd 

b) If ipj = for j odd, we introduce 
(32) 



n —2nikx 



+ Co, < X < 1. 



oo n ^ 

1 (fi 



n+l ^—2'iTikx 



+ Co, < X < 1. 



J even 

Remark. The previous conditions on </?j are not particularly restrictive. 
Indeed, one can easily construct infinitely many polynomial sequences pos- 
sessing the prescribed parity properties. Their generating function has the 
general form 

{e^-l)g{t) ^ n\ 

where g (t) is any real analytic function such that lim(_>o g {t) = 1 and 
satisfying the parity condition g (t) = g (— t) for the case a) and g (t) = 
—g (— t) for the case b). 

Examples 

a) Take 



(e* — 1) sint n\ 

^ ' n=0 



We easily deduce 
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1 1 

ao (x) = 1, ai (x) = x- -, 02 (x) = x'^ - x + -, 

03 a; = X a; H — x , 

^ ^ ^ 2 2 2 

23 

(35) a4(x) -2x^ + 3x2 -2x + —,... 

30 

b) Consider the generating function 
^'^^ (e^-l)cost =2:^-(^);^ 

^ n=0 

We immediately get 

1 7 

6o(x) = l, 6i(x) = x--, 62 (a;) = x^ - X + -, 

2 D 

T/\ 3 *^2 

03 (x = X X H — X , 

^ ^ 2 2 2 

1 7Q 

(37) 64(2;) = x^- 2x3 + 7x2- 6x + —,... 

oU 

4.2. Extended delta rational functions. Let n G Z, g G D. A straight- 
forward generalization of 5n (q) adapted to the chosen Appell polynomials 
is provided by the extended delta rational function A„ (q). 

Definition 11. The extended delta rational function, for any n G Z and 
q €0, is defined by 

Er=i^V n>0 



(38) An iq) = 
where 

n 

(39) ak (n) = ^ 



ET=i^kin)q'' n<0 
1 fi 



7! k'^+^-i ' 

j even or odd 

In eq. (j39p . the summation should be understood either over the even 
values of j or over the odd ones, depending on the choice of the polynomials 
(j3ip or (j32p , respectively. The above definition is motivated by the following 
result, which provides an extension of the construction proposed in Section 
3. Our aim is to obtain the hyper functional equations associated to the 
proposed Appell polynomials. As a consequence of the Fourier expansion 
([30]) and of relations ([3T]) - ([32]) we get the relation between extended delta 
rational functions and Appell sequences. 
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Lemma 12. The following inversion equations, generalizing relation |7^ , 
hold 



(40) A„(g) + (~irA„((7- 

and 
(41) 







if n > 1 



(27ri)-" 



A„ (g) + (-l)"+i A„,(g 



if n>\ 

■^Q-n (^) ^/ n < 



where Pn and Qn are respectively the Appell polynomials \31\) and i33\). to 
which the constant cq has been subtracted. The inversion relations hold for 
all q ^ 1 if n > 0, whereas for n < q can be taken in C\ [0, +oo] . 



By analogy with formulae (jlSp . the associated periodic functions and 
distributions are defined as 



(42) 



Pn{x) 



Pri{x-[x\) 



if n > 1 



-1 + 6j (x) if n = 



jj, (x) if n < 

with an analogous definition for the case of polynomials Qn{x). For all 
n G Z we have the hyperfunctional equations 



(43) 
and 



A„ (,;+) + (-l)"A„, (g 



,-1 



(27ri)-"P_„ (x), 



(44) 



A„ (g+) + (-l)"+i A„ 



_N -1 



(27ri)-"Q„„ {x). 



The proof of these relations is again a direct consequence of the previous 
definitions. Also, by denoting with P„ the restriction to [0, 1] of P„ and 
with if-p the function in (C\ [0, 1]) which represents the hyperfunction 
P„, namely 



(45) 



1 1 



X - T / jo,l] 

we find, for any choice of the sequence {Pn[x)}^^^ and {Qn{x)}^^^ our 
main result, i.e. the generalized Lipschitz formula 



(46) 

Y,^pJr + k) = 2i{27Ti)-'' < 



{ A-n{q) if \q\ < 1, i.e. > 

-1)"-^ A_„ (g-i) if \q\ > 1, i.e. 3t < 
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and the corresponding formula for Q„ (again the usual Lipschitz formula 
corresponds to n < —1). The explicit computation of ip-p^ (r) and LpQ (r) 
is completely analogous to the proposed construction for the Bernoulli poly- 
nomials and we will not repeat it here. 

5. A CONNECTION WITH THE LAZARD FORMAL GROUP 

We describe here briefly an interesting connection between formal groups, 
hyperfunctions, and the so-called universal Bernoulli polynomials. 

Given a commutative ring with identity -R, we will denote by R [xi, X2, .•] 
the ring of formal power series in rri, X2, ... with coefficients in R. Following 
[6] , [1] , we recall that a commutative one-dimensional formal group law over 
R\s a. two-variable formal power series <I> (x, y) £ R [x, y] such that 

1) <^>{x,0) = <^>{0,x) = X 

and 

2) $($(x,y),z) = $(x,$(y,z)). 

When $ (x, y) = ^ {y, x), the formal group is said to be commutative. The 
existence of an inverse formal series (p (x) € R [x] such that $ (x, ip (x)) = 
follows from the previous definition. 

Let us consider the polynomial ring Q[ci,C2,...] and the formal power 
series 

(47) F(s) =S + Ci^+C2y + ... 

Let G (t) be the associated inverse series 

2 3 

(48) G(t)=t-ci| + (3c2-2c2)| + ... 

so that F {G (t)) = t. The series (P7|) and (08]) are called formal group 
logarithm and formal group exponential, respectively. The formal group 
law related to these series is provided by 

$(si,S2) = G(F(si)+F(s2)) 

and it represents the so called hazard's Universal Formal Group It is 
defined over the Lazard ring L, i.e. the subring of Q [ci, C2, ...] generated by 
the coefficients of the power series G {F (si) + F (52))- 

In [IS], the universal Bernoulli polynomials {x,ci, ...,Ck, ■■■) = B^ (x) 
related to the formal group exponential G have been introduced. They are 
defined by 

(49) _i_^e-* = ^ijG(^)^^ 
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The corresponding numbers by construction coincide with the universal 
Bernoulh numbers discovered in [5], and generated by 

(50) o7J)=E«4 -<=^^ 

^ ' n>0 

Observe that when a = I, Ci = then F (s) = log (1 + s) , G (t) = 

e* — 1, and the universal Bernoulli polynomials and numbers reduce to the 
standard ones. Many other examples of Bernoulli-type polynomials consid- 
ered in the literature are obtaining by specializing the rational coefficients 
Cj. The reasons to consider such a generalization of Bernoulli polynomi- 
als are manifold. First, to any choice of the coefficients ci,C2,... it corre- 
sponds a sequence of polynomials of Appell type sharing with the standard 
Bernoulli polynomials many algebraic and combinatorial properties [17]. In 
particular, the associated numbers satisfy universal Clausen-von Staudt and 
Kummer congruences, as shown in the Appendix. Also, in the same way 
as the Riemann zeta function is associated with the Bernoulli polynomials, 
it is possible to relate the polynomials (|49p with a large class of absolutely 
convergent L-series. Their values at negative integers correspond to the 
generalized Bernoulli numbers, and they possess as well several interesting 
number-theoretical properties p^. Observe that the polynomials (j33p be- 
long to the class ()49p . i.e. the coefficient Cj are all rational, if the functions 
g (t) verify the further condition (0) G Q for any n. This is the case, for 
instance, for the proposed examples of the polynomials ([M]) and ([36]) . As 
an immediate consequence, the numbers generated by 

t °° t"^ 

with the prescribed conditions on g (t) do satisfy by construction the Clausen- 
von Staudt and Kummer congruences and many others. These considera- 
tions enable us to associate one-variable formal groups with suitable classes 
of hyper functions, as a consequence of the previous construction. 



6. Appendix: congruences 

The Clausen-von Staudt congruence [8], one of the most beautiful of 
mathematics, states that 

(52) 5„ + V - G Z, 

P 

p\n 

where Bn denotes the n-th Bernoulli number. This proves the strict link 
between Bernoulli numbers and prime numbers. Many generalizations of 
this result have been obtained in the literature in the last decades. In an 
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attempt to clarify the deep connection between these congruences and alge- 
braic topology, in [5] Clarke proposed the notion of universal Bernoulli num- 
bers Bn, defined as (ISOh . and proved the remarkable universal von Staudt 
congruence. 

If n is even, we have 

^ ^n/{p-l) 

(53) Bn = - — mod Z [ci, C2, ...] ; 

p — 1 |n 

p prime 

If n is odd and greater than 1, we have 

(54) Bn = -^^-^ mod Z[ci,C2,...]. 

Like the classical ones, the universal Bernoulli numbers as well play an im- 
portant role in several branches of mathematics, in particular in complex 
cobordism theory (see e.g. [2], and [H]), where the coefficients c„ are iden- 
tified with the cobordism classes of CP". 

Kummer congruences are also relevant in algebraic topology, and in defin- 
ing p-adic extensions of zeta functions. We also have an universal Kummer 
congruence p]. Suppose that n 7^ 0, 1 (mod p — 1). Then 

(55) = — Cp_i mod pZp[ci,C2, ...] . 
n + p — 1 n 

Other related congruences can be found in [1]. 
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